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DYNAMIC STRESS ANALYSIS OF SMOOTH AND NOTCHED FIBER 


COMPOSITE FLEXURAL SPECIMENS 


P. L. N. Murthy and C. C. Chamls 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

A detailed analysis of the dynamic stress field In smooth and notched 
fiber composite (Charpy-type) specimens Is reported In this paper. The analy- 
sis Is performed with the aid of the direct transient response analysis solu- 
tion sequence of MSC/NASTRAN. Three unidirectional composites were chosen for 
the study. They are S-Glass/Epoxy , Kevlar/Epoxy and T-300/Epoxy composite 
systems. The specimens are subjected to an Impact load which Is modeled as a 
triangular Impulse with a maximum of 2000 lb and a duration of 1 ms. The re- 
sults are compared with those of static analysis of the specimens subjected to 
a peak load of 2000 lb. For the geometry and type of materials studied, the 
static analysis results gave close conservative estimates for the dynamic 
stresses. Another Interesting Inference from the study Is that the Impact 
Induced effects are felt by S-Glass/Epoxy specimens sooner than Kevlar/Epoxy 
or T-300/Epoxy specimens. 


INTRODUCTION 

Smooth and notched flexural specimen testing continues to be popular In 
the composites community for characterizing and/or qualifying fiber composites. 
Reasons for the popularity are: (1) simplicity, (2) adaptability to adverse 
environments, (3) availability of simple equations for data reduction, and (4) 
availability of written ASTM standard testing procedures. In addition to 
these reasons, flexural specimen testing forces the material to respond like a 
structure by simultaneously subjecting It to tensile, compressive and shear 
stresses. Furthermore, the notched flexural specimen (Charpy Impact test 
specimen) testing enables easy determination of fracture toughness and Impact 
resistance. 

Fracture In general (be It static, quasi-static or dynamic) Is a dynamic 
event and, as such, Is a very complex process. Fracture progression In a 
flexural specimen is controlled by the local dynamic stress field which Is also 
complex. The local dynamic stress field Is characterized by the stress waves 
which are normal, shear, flexural and surface waves. Each of the dynamic 
stresses can Initiate failure (or a defect) at a point and/or propagate this 
defect to fracture. For a complete understanding of these phenomena a detailed 
description of the local dynamic stress field Is a prerequisite. This requires 
the use of a complex transient analysis Involving direct time Integration as 
opposed to the standard modal synthesis methods such as those available In 
some general purpose finite element structural analysis codes. Solution 27 
(direct transient response) of MSC/NASTRAN Is one of the available tools for 
such analyses. The objective of the present paper Is to report a detailed 
analysis of the dynamic stress field In smooth and notched Charpy specimens. 



Although the notched Charpy test specimen has been used for years In test- 
ing metals and recently In testing composites, no analysis has been performed 
to determine the detailed dynamic stress state variation In the notch vicinity. 
In general, the physical problem of the notched Charpy test specimen Is dynamic 
and non linear; solution of this problem Is difficult. However, a good first 
order approximation may be obtained by assuming linear behavior and a quasi- 
static load. Reference 1 reports such an analysis which uses the static solu- 
tion sequence of COSMIC NASTRAN. The main conclusions of reference 1 are: 1) 

the stress state Is biaxial, and 2) the Charpy test specimen Is not suitable 
for assessing the Impact resistance of nonmetalllc fiber composites directly. 

In the present work the effort Is directed towards understanding the stress 
wave propagation and the attendant dynamic stress field in fiber composite 
smooth and notched flexural specimens. 


ANALYSIS 

In this section, the specimen geometry, the finite element Idealization, 
the finite element analysis method, and the composite systems analyzed are 
described . 


Specimen Geometry 

The geometry of the Charpy test specimen (ASTM STO E23-7) Is shown In 
figure 1. As can be seen In this figure, the overall length of the specimen 
Is 2.164 In. and the length between supports Is 1.574 In. The specimen width 
Is 0.394 In. The specimen unnotched depth Is 0.394 In. and the depth at the 
notch Is 0.315 In. The notch Is 0.079 In. deep and has a 45° opening. 


Finite Element Idealization 

Two finite element Idealizations of Charpy test specimens — one without 
a notch (smooth specimen) and one with a notch were selected for the present 
study. The details of the meshes are shown In figures 2 and 3. The material 
properties are assumed to be uniform, orthotropic, and obey a linear stress- 
strain law throughout the analysis. In addition, the specimen Is assumed to 
be In a state of plane stress. The plane stress assumption Is also justified 
from the physics of the problem. The width restraints at the notch-tip are 
negligible because of the very low value of the respective Poisson's ratio. 

For the present analysis, the plane stress assumption Implies that the stresses 
are permitted to vary along the specimen length and through the thickness but 
not across the width. This reduces the stresses to be calculated to three, 
two normal and one shear. 

With these assumptions, plane stress finite elements can be used to model 
the Charpy test specimen. For the smooth specimen all the elements are quad- 
rilateral. The notched specimen Is modeled with both triangular elements and 
quadrilateral elements. The triangular elements are used as transition ele- 
ments In the areas around the the supports, the load application point and the 
notch. These are the regions where maximum stress concentrations are expected 
to occur and therefore are provided with a finer mesh. The boundary conditions 
prescribed are such that the node at the left support Is constrained from x, y 
and z displacements, and the node at the right support Is constrained from y 
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and z displacements. In addition, for the notched specimen, three nodes clos- 
est to the right support are constrained from displacement In the y direction. 
The specimens are subjected to an Impulse loading. The form of Impact Is a 
triangular function with peak load of 2000 lb occurlng at 500 ys of a total 
contact time of 1000 ys. 

The statistics of the finite element representation for the smooth and 
notched specimens are as follows: 


1) Smooth Specimens. 

Number of nodes or grid points 1147 

Number of displacement degrees of freedom (DOF) (2 degrees 

of freedom per node) 2294 

Number of quadrilateral plate elements (CQUAD4) 1080 

DOF eliminated using the boundary conditions (2 from u = 0 
and v = 0 at the left support and 1 from v=0 at the 

right support) 3 

Number of free DOF (2294 - 3) 2291 

2) Notched Specimens. 

Number of nodes or grid points 656 

Number of displacement degrees of freedom (DOF) 1312 

Number of quadrilateral plate elements (CQUAD4) 544 

Number of triangular plate elements (CTRIA3) 82 

Total number of elements (544 «• 82) 626 

DOF eliminated using the boundary conditions (2 from u = 0 
and v = 0 at the left support and 3 from v = 0 for three 
nodes at the right support) 5 

Number of free DOF (1312 - 5) 1307 


Finite Element Analysis Method 

The MSC/NASTRAN general purpose structural analysis finite element com- 
puter program Is used for the finite element analysis. The specific elements 
used are Identified as CTRIA3 and CQUAD4. They are Isoparametric constant 
strain elements. NASTRAN obtains the solution using a displacement formulation 
via rigid format solution sequence No. 27. This solution sequence employs a 
direct time Integration scheme to obtain the transient response of a structure. 

The solution sequence No. 27 of MSC/NASTRAN uses the Integration algorithm 
based upon the Newmark Beta method (ref. 2). It provides stable results for 
the widest possible spectrum of practical problems without sacrificing either 
accuracy or efficiency. For complete details of the MSC/NASTRAN analysis 
(ref. 3) should be consulted. A brief description Is given In the following 
paragraphs . 

The differential equations of a linear structural problem may be written 
In the general matrix form 


[p 2 M f pB ♦ K] {u} = (F) (1) 

where p = d/dt. M and K are the mass and stiffness matrices, B Is asso- 
ciated with the damping matrix, u Is the vector of displacements and F Is 
the load vector. The numerical Integration Is achieved by replacing p? and 
p by finite difference operators and then using explicit Integration. The 
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outputs of the transient analysis module Include velocities, and accelerations 
as well as displacements. The output can be requested at even multiples of 
the Integration time step. This feature affords some economy In output data 
preparation In cases where small time step Is needed for greater accuracy. 


Composite Systems Analyzed 

Flexural (Charpy-type) test specimens made from three typical composite 
systems are analyzed. They are: T-300/Epoxy, Kevlar/Epoxy, and S-Glass/Epoxy 

composites. The specimens are all unidirectional composites with the fibers 
parallel to the length (x-axIs, fig. 1) of the specimen. 

The plane stress-strain relationship (stiffness) coefficients required to 
Input to NASTRAN are summarized In table I. These properties are obtained by 
using the resident data-bank In the composite mechanics computer code ICAN 
(ref. 4). The relationships between the NASTRAN stiffness coefficients (G's), 
and the usual engineering constants are: 


G 11 

■ W * 1 - 

v ai 2 

u il21 ) 

(2) 

CM 

CD 

= v 9.2i G 11 

= v s.i 2 

6 22 = G 21 

(3) 

G 22 

- - 

w ai2 

W !121 ) 

(4) 

G 33 

= G a.l 2 



(5) 

G 1 3 

= G 23 = G 31 

= G 32 

= 0 

(6) 


The notation In equations (2) to (5) Is as follows: E^n denotes the longi- 

tudinal modulus, E a 22 the transverse modulus, G 112 the shear modulus, w ^-| 2 
the major Poisson's and vj2i the minor Poisson's ratios. For an elastic 
material the two Poisson's ratios are related by the well known relation 

v a.2i = v a.i 2 E a.22 /E a.i l 

RESULTS AND DISCUSSION 
Load conditions 

Two types of load conditions are used In obtaining the results. The first 
Is a static loading where a force of 2000 lb Is applied at the center of the 
specimen on the top surface. The second Is an Impulse loading. This Is 
modeled as a triangular pulse with a peak value of 2000 lb In the middle. The 
pulse Is modeled to last for 1000 pS. The transient response Is, however, 
obtained for three contact time periods (l.e., 3 ms). A separate normal modes 
analysis Is used to determine the first five natural frequencies. 


Natural Frequencies and Normal Modes 

The NASTRAN normal modes analysis module (Solution 3) Is utilized to de- 
termine the first five natural frequencies and the associated mode shapes for 
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the Charpy test specimens. They are shown In figure 4. The results are shown 
In table II for S-Glass/Epoxy specimens. The frequencies are used to determine 
the time periods which aid In determining the Integration time step for the 
transient analysis. A time step of 5 ys Is chosen for the transient analysis 
based upon the time periods shown In table II. 


Displacement and Stress Wave Propagation 


The bulk wave and shear wave velocities are normally much higher compared 
to the flexural wave velocities. In order to capture the characteristics of 
propagation of these waves a much smaller time step of Integration (0.1 ys) Is 
chosen. The output of displacements and stress contours at various time Inter- 
vals are saved and displayed graphically In figures 5 to 10. Figures 5 to 7 
depict the dynamic displacement propagation In S-Glass/Epoxy , Kevlar/Epoxy and 
T-300/Epoxy specimens. Two bulk wave velocity parameters and one shear wave 
velocity parameter are defined below to aid the following discussion of the 
displacement wave propagation results: 


81 1 = 

G n /p 

(8) 

B22 = 


(9) 

SI 2 = 

6 1 2 /p 

(10) 


where p Is the mass density of the material. 

The computed values of Cgn, Cq22> and C S12 are shown In table III for 
the three composite systems under study. The velocities are expressed In 
1n./ys. The trend Indicated by Cb 22 f° r three materials (the transverse 
shock wave travels fastest In S-Glass/Epoxy and slowest In Kevlar/Epoxy) Is 
seen clearly In figures 5 to 7 specifically the frames after 3 and 5 ys. A 
rough estimate of the normal wave velocity can be obtained by counting the 
number of elements that appear to be affected by the Impact from the figures 5 
to 7. The normal wave velocity estimates from the 1 and 3 sec frames are shown 
below: 


Composite system 

Number of 
1 ysec 

elements 
3 ysec 

Velocity 
1 ysec 3 ysec 

Average 

S-Glass/Epoxy 

10 

25 

0.1313 

0.1094 

0.1204 

T-300/Epoxy 

7 

16 

.0788 

.0700 

.0744 

Kevlar/Epoxy 

5 

12 

.0657 

.0520 

.0589 


These values are In close agreement with the theoretical values shown In 
table III under Cg22- The same trend Is also seen for the 3 and 5 y sec 
frames. It can be concluded that the Initial shock travels with the bulk wave 
velocity Cg22 along the direction of Impact. 

Once the normal shock reaches the bottom of the beam, the wavefront ap- 
pears to be moving In the longitudinal direction forming a flexural wave. The 
velocities of the waves traveling In the longitudinal direction can also be 
determined approximately with the same technique mentioned earlier. The fol- 
lowing are the details for the frame after 13 ys: 
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Composite system Number of 

elements 


Velocity 


S-Glass/Epoxy 11 0.0510 In/ysec 

T-300/Epoxy 9 .0416 In/ysec 

Kevlar/Epoxy 7 .0320 In/ysec 

The above velocities appear to have the same trend as depicted by the wave 
velocity parameter Csi 2 shown table III. However, the waves along the 
longitudinal direction appear to move significantly slower than that Indicated 
by Csi2- This Is probably due to the coupling between the flexural wave and 
the shear wave. The flexural wave velocity Is significantly slower than the 
shear wave velocity. For example, the smooth S-Glass/Epoxy specimen under 
study has a flexural wave velocity given by 

C p = 2fi.f = 0.0137 In/ysec 

where 8. Is the length between the supports, f Is the first fundamental 

frequency. (It Is assumed that the wave number Is 1 and the beam deflects Into 

a half wave. ) 

The transient stress response after one percent (10 ys) of the contact 
time are shown In figures 8 to 10 In the form of stress contours. All stresses 
are localized at this early time. The longitudinal stress (on) In the 
Kevlar/Epoxy and T-300/Epoxy composites Is of about the same magnitude while 
that In the S-Glass/Epoxy Is about half as much. The normal stress ( 022 ) and 
shear stress (o-j 2 ) are of about the same magnitude for all three composite sys- 
tems. It appears from these stress results that, under the same Impact condi- 

tions, the stress In the Kevlar/Epoxy and the T-300/Epoxy will reach fiber 
fracture stress magnitudes considerably earlier (about half the time) than In 
the S-Glass/Epoxy composite. Two Implications follow relative to the same 
stress magnitude: (1) the S-Glass/Epoxy composite will sustain greater Impact 
load prior to fracture than the T-300/Epoxy, and (2) the rapid compressive 
stress built up will cause longitudinal compression failure accompanied by 
substantial bending deflection In the Kevlar/Epoxy composite thus Increasing 
the Impact required to Induce fracture. 


Transient Response of Smooth Specimens 

The displacements, velocities, accelerations and stresses are obtained for 
a total contact time of 1000 ys for S-Glass/Epoxy smooth specimens. In these 
computations a time step of 5 ys Is used. The output Is saved for every two 
time steps. Two points A and B as shown In figures 2 and 3 are selected for 
study. A Is the load point and B Is the opposite point at the bottom for 
smooth and at the notch-tip for notched specimens. 

The transient response results appear In figures 11 to 16. Figures 11 to 
13 show the longitudinal (axial) and the transverse (bending/flexural) com- 
ponents of displacement, velocity and acceleration of point A plotted against 
time. The corresponding results for point B are not shown as they are similar 
to point A response both qualitatively and quantitatively. It can be concluded 
from these figures that the response Is primarily In the first flexural mode. 
For example (from fig. 11) the number of cycles In 1 ms Is counted as a little 
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over 6. From table II the first fundamental time period Is 161.3 ys which 
Implies that 6.2 (1000/161.3) cycles of response In first mode can be expected. 

The stress response results are shown graphically In figures 14 to 16. 

Each figure contains the stress response of elements near points A and B. The 
longitudinal stress behavior appear to be primarily the first flexural mode 
response. However, the magnitude of stress near point A Is about one and a 
half times greater than the magnitude of stress near point B. This Is to be 
expected because of stress concentration near the load point A. The transverse 
normal ( 022 ) and the shear ( 012 ) stress Increase linearly to a peak when the 
load takes the maximum value of 2000 lb and then decrease linearly to a zero. 
The corresponding stresses for the element near point B are Insignificant and 
therefore can not be shown distinctly when drawn to the same scale. 

The response after the load removal Is not shown In the above figures. 

The computations, however, are conducted for three contact time periods. It 
Is found that, for all practical purposes, the specimen remains static. A 
typical response curve Is shown In figure 17 for point A. 


Static and Dynamic Stress Contours 

Smooth specimen . - The static stress contours under a load of 2000 lb and 
the dynamic stress contours at the peak load of 2000 lb are shown In figures 
18 to 20. They appear to be Identical. The peak values for the dynamic load 
case, however, are slightly lower than those for the static case. While the 
maxima for the longitudinal stress ( cr-j -j ) occur at the center of top and the 
bottom, the transverse normal stress and the transverse shear stress have the 
maxima almost adjacent to each other near the load application point. The 
transverse shear stress and the axial stress distribution Is seen to approach 
the classical Euler-Bernoul 1 1 theory predictions away from the load application 
point. Steep stress gradients are observed near the load application point. 

Notched specimen . - The results for notched S-Glass/Epoxy specimen are 
shown In figures 21 to 23. The static and dynamic stress contours at the peak 
load show similar trend as observed In the case of smooth specimens. Figure 24 
shows the stress Intensity near notch-tip In relation to the far field 
stresses. This distance from the notch-tip Is measured along the longitudinal 
direction towards the left support. It Is seen that all the three stresses 
attain very high peaks, Indicating severe local stress Intensities near the 
notch-tip. As the stress allowables for 022 and <*12 are generally an order 
of magnitude lower than that for oyj, one can expect a matrix Initiated fail- 
ure followed by fiber fractures at this location. 


GENERAL DISCUSSION 

The response analysis of smooth and notched Charpy-type flexural specimens 
subjected to a triangular Impulsive load provides Insight Into the nature of 
stress wave propagation, failure mechanisms and the relation between the static 
and dynamic responses. The contact time of the Impact load Is approximately 6 
time periods of the first flexural mode. The bulk wave velocities Cgyi and 
C B22* and the sf)ear wave velocity Csi 2 are very high compared to the flexural 
wave velocity for the materials and the geometry under study. Hence, the re- 
sponse Is observed at early times of the Impact event (of the order of a micro- 
second). It shows that the Initial waves travel with a bulk wave velocity 
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C B22 till they reach the bottom surface of the specimen. The waves then 
appear to travel longitudinally towards the supports. The approximate calcu- 
lations based upon the deformations at various times Indicated that the wave 
velocity Is significantly lower than the shear wave velocity but much 

higher than the flexural wave velocity. The wavefront appears to Induce a 
flexural wave as It progresses towards the supports. 

At longer times, the response Is primarily In the first flexural mode. A 
comparison with the static response Indicates that little or no difference 
exists In the magnitude for stresses at the peak load. The static predictions, 
however, are on the conservative side. 

The response of notched fiber composite specimens show severe stress In- 
tensities near the notch region Indicating local failures In shear and trans- 
verse tension Initially followed by fiber fractures since the longitudinal 
stresses approach fracture stress magnitudes. 

The response after the load removal reduces to mere noise; the structure 
remains practically static. 


SUMMARY OF RESULTS 

The results of the study on static, dynamic and transient response of 
smooth and notched, fiber composite, Charpy specimens are listed below: 

(1) The wave propagation velocities can be estimated from the early time 
displacement propagation response; the estimates are In fair agreement with 
the theoretical predictions. 

(2) The first five flexural mode shapes and frequencies show coupling 
effects from the thickness stretch and the transverse shear type modes; the 
frequencies are not Integer multiples of the fundamental frequency. 

(3) The dynamic and static peak load stress contours are almost Identi- 
cal. The static peak load stress magnitudes are slightly higher. Stress pre- 
dictions based upon a quasi-static approach lead to conservative estimates. 

(4) The transient displacement and velocity response appear to be pri- 
marily In the first flexural mode. The acceleration response shows contribu- 
tions from higher modes. 

(5) The structure responds for the load duration time period only. The 
steady state response after the load removal appears to be negligible. 

(6) The transverse normal and shear stresses vary linearly with time and 
follow the load path. 

(7) The notch-tip region develops severe stress concentrations and any 
of the three stresses could cause or Initiate a failure. 

(8) Based upon the transient stress response, It appears that the fail- 
ures In Charpy specimens are Initiated at the notch-tip by the shear stresses 
and the transverse tensile stresses followed by fiber fractures. However, the 
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failures In smooth specimens are probably Initiated by a combined stress state 
near the load application point followed by local Interply delamlnatlons . 
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TABLE I. - NASTRAN PLANE STRESS - STRAIN RELATIONSHIPS 


Composite 

Mass density 

Stress-strain coefficient 

Orthotropy Ratios 

system 

( xlO 4 ) 

( lb. sec 2 /1n 4 ) 

G 11 

G 12 “ G 21 

G 22 

(xlO 6 psl ) 

G 33 

G ll /G 22 

G ll /G 33 

S-Glass 

/Epoxy 

1 .813 

8.7900 

1.1700 

3.2600 

1.2700 

2.70 

6.92 

Kevlar 

/Epoxy 

1.249 

12.385 

.20082 

.57377 

.31635 

21.59 

39.15 

T-300 

/Epoxy 

1.405 

17.896 

.30275 

1.1318 

.50871 

15.81 

35.18 


TABLE II. - NATURAL FREQUENCIES AND TIME 
PERIODS FOR S-GLASS/EPOXY SMOOTH AND 
NOTCHED SPECIMENS 


Specimen 

type 

Mode 

number 

- 

Frequency 

(cycles/sec) 

Time period 
(ps) 

Smooth 

1 

6200.67 

161.3 



11659.54 

85.8 



19993.02 

50.2 



31062.8 

32.2 



43913.02 

22.8 

Notched 

1 

9030.16 

110.7 



14631 .33 

68.3 



24596.56 

40.7 



27775.86 

36.0 


mm 

36685.27 

27.3 

























TABLE III. - WAVE VELOCITY PARAMETERS 
FOR THE SELECTED COMPOSITE SYSTEMS 
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Figure 1. - Geometry of ASTM Charpy test specimen. 
(All dimensions in inches. ) 


P(t) 



Figure 2. - Finite element idealization of smooth specimen. 



Figure 3. - Finite element idealization of notched Charpy 
specimen. 



L f ■ 6200. 67 Hz. 



1. f = 9030. 157 Hz. 



3. f = 19993. 02 Hz. 



4. f • 31062. 80 Hz. 



4. f - 27775. 86 Hz 



5. f « 43913. 02 Hz. 5. f = 36685. 27 Hz. 

Figure 4. - Natural frequencies and mode shapes for 
S-glass/epoxy smooth and notched flexural specimens. 
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Figure 5. - Dynamic displacement propagation in S-glass / 
epoxy composite flexural specimen. (Triangular impulse: 
2000 lb maximum 1000 nsec duration. ) 
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Figured - Dynamic displacement propagation in Kevlar/ 
epoxy composite flexural specimen. (Triangular impulse: 
2000 lb maximum 1000 nsec duration. ) 
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Figure 7. - Dynamic displacement propagation in T-300/ 
epoxy composite flexural specimen. (Triangular impulse: 
2000 lb maximum 1000 psec duration. ) 


SYMBOL VALUES, psi 

1 -L 328899E+03 

2 -L 143674E+03 

3 -9. 584489E+02 

4 -7. 732237E+02 

5 - 5, 879985E+02 

6 -4 027733E+02 

7 -2. 173482E+02 

8 -3. 232298E+01 

9 1. 529022E+02 

10 3. 381274E+02 



(a) ojj Contours. 


1 -L052444E+03 

2 -9.068754E+02 

3 -7. 613066E+02 

4 -6. 157379E+02 

5 -4701691E+02 

6 3. 246003E+02 

7 1. 790316E+Q2 

8 3. 346283E+01 

9 1. 121059E+02 



(b) 022 Contours. 


1 -5. 760680E+02 

2 -4 970201E+02 

3 -4 120802E+02 

4 -3.689114E+02 

5 -2. 998825E+02 

6 -1. 808336E+02 

7 -1. 017847E+02 

8 -2. 273584E+01 

9 5. 631304E+01 
10 1. 353619E+02 



(c) Oj 2 Contours. 


Figure 8. - Dynamic stress contours in 
S -glass/epoxy composite flexural 
specimen after 1 percent contact time. 


SYMBOL VALUES, psi 

1 -2. 362948E+03 

2 -2. 065930E+03 

3 -L768912E+03 

4 -L471895E+03 

5 -1. 174877E+03 

6 -8 778589E+02 

7 -5. 308411E+02 

8 -2 838233E+02 

9 1 319453E+01 

10 3. 102123E+02 



(a) On Contours. 


1 -1. 913533E+03 

2 -a 775216E+02 

3 -7. 415100E+02 

4 -6. 054964E+02 

5 -4 694862E+02 

6 -3. 334752E+02 

7 -L 974635E+02 

8 -6. 145193E+01 

9 7. 455969E+01 



(b) o 2 2 Contours. 


1 7. 771991E+02 

2 -6. 946070E+02 

3 - 5. 920160E+02 

4 -4 994229E+02 

5 -4 068309E+02 

6 -3. 142388E+02 

7 -2. 216468E+02 

8 -1. 290548E+02 

9 -3. 646273E+01 

10 5. 612931E+01 



(c) Oj 2 Contours. 


Figure 9. - Dynamic stress contours in 
Kevlar/epoxy composite flexural specimen 
after 1 percent of the contact time. 
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Figure 10 - Dynamic stress contours in 
T-300/epoxy composite flexural specimen 
after 1 percent contact time. 
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Figure 12. - Velocity response of 
point A for S-glass/epoxy smooth 
flexural specimen. 
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Figure 13. - Acceleration response 
for S-glass/epoxy smooth flexural 
specimen. 



Figure 14 - Axial stress o xx re- 
sponse of top and bottom elements 
near points A and B for S -glass/ 
epoxy smooth flexural specimen. 



Figure 15. - Transverse normal 
stress Oyy response of top and 
bottom elements near points A 
and B for S-glass/epoxy smooth 
flexural specimen. 



Figure 16. - Transverse shear stress 
o xy response of top and bottom 
elements near points A and B for 
S-glass/epoxy smooth flexural 
specimen. 




Figure 17. - Displacement response 
of point A for three contact time 
periods of S-glass/epoxy smooth 
flexural specimen. 
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(a) Static stress contours. 
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(b) Dynamic stress contours. 


Figure 1& - Comparison of static and dynamic 
axial stress o xx contours at peak load for 
S -glass /epoxy smooth flexural specimens. 


SYMBOL VALUES, psi 

1 -5. 259623E+04 

2 -4. 533239E+04 

3 -3. 806855E+04 

4 -3. 080471E+04 

5 -2. 354087E+04 

6 -1. 627703E+04 

7 -9. 013191E+03 

8 -1. 749351E+03 

9 5. 514490E+03 

10 1. 277833E+04 



la) Static stress contours. 
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(b) Dynamic stress contours. 


Figure 19. - Comparison of static and dynamic 
transverse normal stress o„y contours at 
peak load for S -glass /epoxy smooth flexural 
specimen. 
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Figure 20. - Comparison of static and dynamic 
transverse shear stress o xy contours at 
peak load for S-glass/epoxy smooth flexural 
specimen. 
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(a) Static stress contours. 
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(b) Dynamic stress contours. 


Figure 21. - Comparison of static and 
dynamic axial stress o xx contours 
near the notch region at peak load 
for S-glass/epoxy notched flexural 
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(a) Static stress contours. 
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(b) Dynamic stress contours. 


Figure 22. - Comparison of static and 
dynamic transverse normal stress 
Oyy contours near the notch region 
ar peak load for S-glass/epoxy 
notched flexural specimen. 
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(a) Static stress contours. 
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(b) Dynamic stress contours. 


Figure 23. - Comparison of static and 
dynamic transverse shear stress 
contours near the notch region 
at peak load for S-glass/epoxy 
notched flexural specimen. 
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